High-performance diamond-based single-photon sources for quantum communication 
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Quantum communication places stringent requirements on single-photon sources. Here we report a 
theoretical study of the cavity Purcell enhancement of two diamond point defects, the nickel- nitrogen 
(NE8) and silicon-vacancy (SiV) centers, for high-performance, near on-demand single-photon gen- 
eration. By coupling the centers strongly to high- finesse optical photonic-bandgap cavities with 
modest quality factor Q — O(IO^) and small mode volume V — O(A^), these system can deliver pi- 
cosecond single-photon pulses at their zero-phonon lines with probabilities of 0.954 (NE8) and 0.812 
(SiV) under a realistic optical excitation scheme. The undesirable blinking effect due to transitions 
via metastable states can also be suppressed with O(10~^) bUnking probability. We analyze the 
application of these enhanced centers, including the previously-studied cavity-enhanced nitrogen- 
vacancy (NV) center, to long-distance BB84 quantum key distribution (QKD) in fiber-based, open- 
air terrestrial and satellite- ground setups. In this comparative study, we show that they can deliver 
performance comparable with decoy state implementation with weak coherent sources, and are most 
suitable for open-air communication. 

PACS numbers: 03.67.Hk, 42.50.Dv, 42.60.Da, 76.30.Mi 



I. INTRODUCTION 



Light sources that reliably prepare single photons in 
well-defined spatiotemporal modes on demand are en- 
ablers for experiments in quantum optics py], precision 
optical measurement [2|, [3[ and numerous prospective 
quantum information technologies such as quantum tele- 
portation [4], quantum repeaters [H, and linear-optical 
quantum computers Q. Their most immediate applica- 
tion is in quantum cryptography Q. The faithful and 
practical implementation of the well-known BB84 quan- 
tum key distribution (QKD) scheme [8] relies on the use 
of a high-performance true single-photon source (SPS). 
Most QKD systems realized to date employ weak co- 
herent sources (WCS) with relatively low single-photon 
throughput. While the resultant key creation rate can 
be improved drastically with decoy states [TO, [H!] , fur- 
ther enhancement and simplification of practical BB84 
implementations are possible using SPSs [12]. 

There are numerous schemes for single-photon gen- 
eration including two-photon down-conversion 13|, 
Coulomb blockade of electrons state reduction lEj , 
and fluorescence of a wide range of quantum systems, 
e.g. atoms [16], trapped ions [17], molecules [18, 19], 
and quantum dots |20|. Various optically-active point 
defects in diamonds are also superb SPS candidates and 
have been extensively investigated. In particular, anti- 
bunching experiments with isolated centers, namely the 
negatively-charged nitrogen- vacancy (NV) [21], nickel- 
nitrogen (NE8) [23, m and silicon-vacancy (SiV) H 
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centers, have also found them to be photostable at room 
temperature and robust against photobleaching. Clock- 
trigger^ NV and NE8 based SPSs were demonstrated in 
Refs. j25|, l26| (for QKD) and [23] respectively. Coupling 
of a diamond nanocrystal contai ning NV to an optical 
fiber has also been reported [23, [l^ and fiber-coupled 
sources are now commercially available [29| . 

Despite these successes, the bare centers are not ideal 
for practical QKD and prospective technologies. Their 
emission is nondirectional and this impedes collection ef- 
ficiency. The nanosecond photoluminescence lifetimes of 
these centers limit their theoretical operational speeds to 
< 500 MHz. The presence of a metastable state leads 
to undesirable blinking and > 100 ns long dark intervals, 
and decays via upper-phonon transitions lead to broad- 
ened emission bandwidths [22, i2J, [30|]- Also, the SiV 
center suffers from fast nonradiative decay that degrades 
its fluorescence [31] and ultimately reliability. Although 
there are other newer, uncharacterized defects [32l, [sl, [sj] 
with large transition dipole moments, we focus on these 
most commonly discussed centers. 

To realize a directional and high-performance source, 
the emitter can be placed inside an optical resonator. 
In this cavity quantum electrodynamic (QED) arrange- 
ment, coupling between an NV center and an active cav- 
ity mode has been studied previously in Ref. [35[ , showing 
dramatic improvements in the spectral purity and emis- 
sion rate. In the first part of this work (Sec. II), we 
study the use of cavities to enhance the emission of the 
NE8 and SiV centers - but with a clear distinction owing 
to appreciable differences in their electronic and vibronic 
structures. We show that wavelength-sized cavities with 
quality factor Q = O(IO^) improve their emission rates 
by over two orders of magnitude, and overcome the pho- 
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toluminescence problem of the SiV center and the com- 
mon bhnking effect. In Sec. Ill, we explore the potential 
of these sources (including the NV) toward various QKD 
implementations. Among silica fiber-based, open-air ter- 
restrial and satellite-ground arrangements, we find that 
the cavity-enhanced centers deliver performance compa- 
rable with the decoy state, and are most suitable for 
open-air communication. 

Photonic crystals (PCs) are ideally suited for realiz- 
ing the required cavity properties. The schematic in 
Fig. [Tfb) envisages such an implementation. A defect in 
a planar photonic-bandgap (PBG) structure defines an 
excellent cavity, and extremely high-Q/l^ cavities [mode 
volume V = O(A^), Q = O(IO^)] for near-infrared (near- 
IR) operation have been demonstrated in silicon [36j]. 
For enhancing the visible emission from these centers, 
compatible hi^h-Q/V cavity designs are also available 
in diamond [33, [^, [s^, HO, [4l| and gallium phosphide 
(GaP) [42|. Although fabrication of cavity-enhanced cen- 
ters in integrated photonics [10] remains challenging, ex- 
perimental demonstration of diamond-based PBG cavity 
of Q = 585, V = 0{X^) pi] and coupling an NV to a 
silicon nitride cavity of Q = 3200 |45|, give cause for op- 
timism. Realizing a coupled NV-cavity system has also 
been achieved with Q > 25000 microdisk in GaP [10]. 
Finally, we note that novel slot- waveguide geometries, 
combined with mirrors, PBG or distributed Bragg re- 
flectors in a Fabry-Perot arrangement, are promising for 
achieving ultra-high optical confinement [43, S, SO] • 



II. CAVITY-ENHANCED DIAMOND CENTERS 

Diamond color centers are vibronic systems that con- 
sist of a ground state manifold of vibrational sublevels 
{\gi)} and an excited state manifold of sublevels {|ej)}, 
where z,j = 0,1,2,... label these levels with increasing 
energy. Transitions between the manifolds are electric 
dipole in nature, giving rise to multiple phonon lines, de- 
noted iPL, in the emission spectrum: |eo) \go) transi- 
tion is the sharp zero-phonon line (ZPL), and the vibronic 
sidebands due to phonon- assisted transitions |eo) \gk) 
(k > 0) are broadened by rapid phononic relaxations. Be- 
low, we rename |eo) as |e) as shorthand. The NE8 center 
has only very weak sidebands at room temperature and a 
prominent ZPL centered at Azpl = 794 nm. The relative 
integrated intensity of the ZPL to the entire spectrum is 
0.7 [22|. For the SiV center, its ZPL centered at 738 nm 
makes up ~ 80% of the emission intensity [50]. Their 
IPLs are roughly centered at 20-30 nm from the ZPLs. 
Thus, we adopt a common three-state model {|^o), 
|e)} with photoluminescence lifetimes r of 11.5 ns (NE8) 
and 2.7 ns (SiV). Added to our model is a nonradiative 
metastable state \h) that is responsible for dark intervals 
(blinking) and the photon-bunching signature in the au- 
tocorrelation function ^^^^ of their emissions [22|,[23]. The 
complete model inside a one-sided, single-moded cavity 
is depicted in Fig. [H 




FIG. 1: (Color online) (a) Schematic of a cavity-center sys- 
tem for single-photon generation. The NE8 and SiV centers 
are modelled as a vibronic system that consists of a single 
excited state |e) and a ground state with two vibrational 
sublevels \gi). \h) is an additional metastable level. The 
center is pumped with an external classical field r(t) acting 
as the trigger pulse. The |e) — \gi) transition is coupled to 
a one-sided, single-moded cavity with coupling strength Qi. 
Dashed arrows denote radiative decay and dotted arrows non- 
radiative relaxations, (b) Proposed design in diamond-based 
PBG lattice, missing holes constitute the cavities and waveg- 
uide. Solid circle (in box) denotes a center excited optically 
by a green laser, (c) NE8 complex in the carbon lattice, con- 
sisting of a Ni atom surrounded by four N atoms, (d) SiV 
structure consisting a substitutional Si atom and an adjacent 
vacancy. 

In dipole and rotating- wave approximations, the 
Hamiltonian of the joint center-cavity system is 

H/h = ^U;aCraa^^cCilcic^ ^ l^i (ajcr^^e + h-C) , (1) 
a i=0,l 

where aaf3 = is the atomic projection operator, \a) 

and \P) denote the atomic states, and al (ac) is the cre- 
ation (annihilation) operator for the cavity mode-c. hcOa 
is the energy of the atomic state \a) and hcOc = 27Thc/Xc is 
the resonant energy of the cavity. The dipole transitions 
iFL are coupled to the cavity mode with single-photon 
Rabi frequency Qi = [^c/(2^eol^)]^/^, parameterized 
by the transition dipole moment di and the cavity mode 
volume V. 

The dynamics of the system coupled to an adjacent 
waveguide (mode-it;) satisfies the master equation for its 
density operator p, 

i 



3 



i=0,l 

where the Liouvihean, defined for an operator A, 

C[A] = ApA^ - 1 (^A^Ap + pA^a) , (3) 

is used to describe the relevant atomic relaxation pro- 
cesses, namely the relaxation to the metastable state 
of lifetime 1/jh with rate 7e, the phononic relaxation 
within the manifold with rate 7^ ^ 7i, the atomic (ra- 
diative) spontaneous transitions iPL with decay rate 7^, 
and other nonradiative decay processes with combined 
rate 7. = uJc/{2Q) is the rate of cavity loss into 
mode-iL' associated with the creation operator aj^. We 
assume that this effect dominates the determination of 
Q. Finally, we include a pulsed excitation with pumping 
rate r{t) that serves as a trigger for single-photon emis- 
sion. The magnitude of the transition dipole moment di 
of iFL transition can be determined using the relation 
(if = STrheoc^ / (jincof) ^ where n ^ 2.4 is the refractive 
index of diamond and hcoi is transition energy. Allow- 
ing the possibility of multiple excitations/de-excitations 
per pulse, we adopt a basis set {|e), l^o), Is'i), |^)} ^ 
{|Oe),|lc),...,|A^c)} ^ {|0^),|l^),...,|iV^)}, where |n,) 
(In^)) denotes the number state of mode-c {w). 

Although the shelving level reduces photon emission, 
the quantum yield of the bare NE8 center remains near 
unity, i.e. 7 = 0. In comparison, the SiV center has 
a very low quantum yield of ~ 0.05 at room tempera- 
ture [5l| largely due to nonradiative transitions so that 
7 = 19/r. Since the branching ratio 70/71 is related to 
the relative integrated intensity of the ZPL to the to- 
tal emission, we take 70/71 = 7/3 (NE8) and 4 (SiV). 
Fitting to the measured autocorrelation function ^^^^ in 
Ref. [22] gives 7e/27r = 17 MHz and jh/^tt = 6.1 MHz 
for the NE8. For the SiV, 7e/27r - 30 MHz and 
7^/27r ~ 10 MHz [2^. In this model, we have explic- 
itly ignored thermal broadening but it can be introduced 
phenomenologically for a more realistic estimate of the 
spectral width of the output wave packet. 

Reliable single-photon generation from the cavity- 
center system requires that losses via atomic relaxation 
and the shelving state be minimized, thus the atomic ex- 
citation should be outcoupled into mode-w faster than 
these processes. In the strong Purcell regime [52] where 
the cavity is resonant with the Mth phonon transition 
and hz > Qm ^ {7i(Vi), 7e, 7/1, 7}, the Rabi frequency 
dominates the rates of the incoherent processes so that 
the state vector evolves as \e,Oc^^w) I^M^lcO^^) 
I^M, Oc, Iw). Notably, tz > ensures rapid outcoupling 
from the cavity mode and suppresses the vacuum Rabi 
oscillations, which otherwise lead to unwanted Rabi rem- 
nants - "ringing" spectral features at the cavity output. 

Here we focus on enhancing the ZPL transition with 
diamond-based PEG cavities of V = {Xc/n)^, which is 
of the order typically obtained in proposed cavity de- 
signs [41]. Although the theoretical limit for the funda- 
mental mode of such a cavity is (Ac/2n)^, these modes 



are normally not used because they are lossy [36|. For 
the upper phonon transitions to be enhanced, a large cou- 
pling Qm ^ Jg and hence sub-wavelength confinement 
[V <C (Xc/n)^] are required due to their associated broad 
linewidth [35[. 

Before presenting the simulation results, the optimal 
operating regime warrants further discussion. The cavity 
enhancement of a quantum emitter is measured by the 
Purcell factor Fp, defined as the ratio of the emission rate 
to the cavity mode to the no-cavity case [53|, [13, [EH . In 
this case, it is given by 

r < Q (4) 

^ 7total^ ^eo7total V 

where 7totai = 27r/r + 7 + 7e is the total decay rate from 
state |e) expressed in units of angular frequency. Purcell 
enhancement is maximized with a high-Q/l^ cavity (i.e. 
small K,)^ but an arbitrarily large Q delays cavity output. 
For a passive structure, there is an upper bound for Q, 
or equivalently, hz > 21^0 • Hence the Q factor is optimal 
when hz ^ 2.51^0 and the Rabi remnants suppressed. 

When the center is pumped too strongly, the use of 
cavities makes multi-photon emission events per trigger 
possible. The excitation parameter r{t) must be op- 
timized to maximize excitation probability and single- 
photon emission probability Pi while suppressing multi- 
photon probability Pm- In general, the temporal extent 
of an excitation pulse must be much shorter than the 
time scales for Rabi oscillation and cavity decay. During 
excitation we can then ignore the effect of decoherence 
and treat the centers are two-level systems. Assuming 
a top-hat trigger pulse of length T and absorption rate 
r, we solve for state populations in the zero, one and 
multi-excitation manifolds. As decoherence only reduces 
populations in finite excitation manifolds, these results 
become the upper bounds for Pi and Pm given by 

-wny^^^-ir^-wnl)], (5) 

Pm < l-Pi-exp(-rT). (6) 

In general, the optimal regime is accessible by using 
ultrashort and intense pulses, but such optical pumping 
may effect changes in the centers. In a femtosecond il- 
lumination experiment by Dumeige et a/., they showed 
photochromism of a NV under a peak intensity of near 
50 GW/cm^ with 55 pJ, 150 fs pulses [56], and no damage 
to the diamond was reported. As similar studies of the 
NE8 and SiV are not present, we assume that these cen- 
ters remain stable for a peak intensity / < 50 GW/cm^. 

Following Eqs. [5] and [6l we choose parameters (T, r) 
= (0.16 ps, 20 THz) and (0.08 ps, 40 THz) for the 
NE8 and SiV respectively. These values correspond to 
/ = Ti? ce^r'^ / (2(f) ^ 40 GW/cm^ and the energy per 
pulse is < 50 pJ. Consequently, emissions from the NE8 
and SiV follow Pi < 0.97 and 0.95 respectively, with 
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Pm = O(10~^). On the other hand, for apphcations 
where minimizing is favored over near unit Pi (e.g., 
QKD in Sec. Ill), shorter pulses can be used. For in- 
stance, for both centers, 20 — 40 fs pulses are suitable for 
achieving Pm = O(10~^) for Pi ~ 0.56. Development of 
modern mode-locking techniques made possible the gen- 
eration of such ultrashort pulses (see Refs. [s?, ^sJIHO] 
for review articles on femtosecond lasers) with high repe- 
tition rates and notably, a titanium-sapphire laser capa- 
ble of emitting 42 fs pulses at a 10 GHz repetition rate 
has been demonstrated [6l|. Gigahertz repetition rates 
can also be achieved by gating multiple ultrashort lasers. 

Armed with this understanding, we present numeri- 
cal simulations of a single excitation cycle. Figure [21(a) 
plots the temporal lapse and envelope of the outcoupled 
photon pu,w = (S'cOc, lw|p|^o,Oc, 1^) in response to an 
excitation pulse that yields Pm = O(10~^). Excitations 
with shorter pulses that yield Pm = 0(10"'') are also 
considered. The cavity-Q values are modest: 3700 (NE8) 
and 1800 (SiV), and they are near optimal as shown in 
Fig.[2^b) where simulations are performed for different Q. 
In the weak cavity regime of > 2Qo^ Purcell enhance- 
ment diminishes with Q. In the strong cavity regime 
where Q exceeds ujc/^^o^ the emission yield via the cav- 
ity channel reduces as the time scale for cavity relaxation 
becomes much longer than the other decay processes. 

The ratios Q/V used in Fig. [2](a) imply emission en- 
hancements with factors P^ = 311 (NEB) and 10 (SiV), 
and Pi = Fp/{1 + Fp) = 0.997 and 0.908 respec- 
tively if excitations are ideal. In the full treatment, 
we expect some reductions - at large times the prob- 
ability of emitting a photon of ZPL, given by Pi = 
{go, Oc, 1^ \p\go, Oc, 1^), is 0.954 for the NEB. For the SiV, 
the effect of fast nonradiative decay 7 in the SiV is still 
significant but the cavity has improved Pi from 0.05 of 
a bare SiV to 0.812. When shorter excitation pulses are 
used. Pi = 0.565 (NFS) and 0.470 (SiV). Irrespective of 
the pulse length T, the mean times for photon output 
are predicted to be 10 ps (NE8) and 5 ps (SiV), nearly 
three orders of magnitude faster than their nanosecond 
free-space lifetimes. 

The cavity-enhanced centers are now capable of op- 
erating at repetition rates of = 30 GHz (NE8) and 
65 GHz (SiV). In these systems, multi-photon events are 
4-6 orders of magnitude less likely than Pi, the proba- 
bility of a dark interval per cycle is O(10~^), and that 
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of emitting a IPL photon is (^1, Oc, lw|p|^i, 



0(10-^) (NE8) and 0(10"^) (SiV). Fourier transforms 
of the temporal profiles yield emission spectrums cen- 
tered at the ZPL with effective linewidth of 0.15 nm 
(NE8) and 0.3 nm (SiV) at zero temperature. Although 
the pulse is non-Gaussian, techniques such as Stark tun- 
ing and Q-switching can be used to optimize pulse pro- 
files [oil, Isl, m, ^ and improve photon indistinguisha- 
bility. Next the full model is solved for different excita- 
tion parameters in Figs. [2l(c,d). The results are in agree- 
ment with Eq. [5]and show suitable (T, r) range for SPS. 
Although small improvements to the emission probabili- 
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FIG. 2: (Color online) (a) Cavity-enhanced photon emission 
at ZPL wavelengths from two diamond optical centers: NFS 
(Ae = 794 nm) and SiV (738 nm). We take V = {Xc/nf and 
Q = 3700 (NFS), 1800 (SiV). The solid curves are calculated 
using top-hat excitation parameters (T, r) = (0.16 ps, 20 THz) 
and (0.08 ps, 40 THz) for the respective centers. The dashed 
curves use shorter excitation pulses of T = 40 fs (NF8) and 
20 fs (SiV). Single-photon emission probability Pi as a func- 
tion of (b) Q, and (c,d) the excitation parameters. Circles 
(squares) indicate the parameters that effect Pm = O(10~^) 
[0(10-^)]. 
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FIG. 3: Photon autocorrelation histogram of emissions from 
the cavity-center system implemented with an (a) NF8 and 
(b) a SiV center. Using quantum trajectory approach, HBT 
experiment is simulated using a pulsed excitation, correspond- 
ing to a top-hat pulse at respective repetition rates 10 GHz 
and 20 GHz for over 5000 emission cycles. 



ties, Pi = 0.995 (NE8) and 0.854 (SiV) while maintain- 
ing Pm = O(10~^), are possible with T ~ 10 fs and 
r ~ 1 PHz pulses, the corresponding intensity is greater 
than those for observing photochromism. 

Photon autocorrelation obtained in a Hanbury-Brown- 
Twiss (HBT) experiment is a test for single-photon 
emission. We simulate the experiment over 5000 
excitation-emission cycles using the quantum trajectory 
approach [66, _67.] . Assuming the use of photon detection 
units with time bin sizes of 0.6 ps (for NE8) and 0.1 ps 
(SiV), the coincidence histograms cx g^'^\r) in Fig. [3] are 
compatible with perfect photon antibunching - suppres- 
sion of coincidence counts (i.e. multi-photon events) at 
and around zero delay. We summarize the performance 
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2.95 


64000 
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120 


0.95 
0.74 
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(•) 
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2.37 
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9 
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1 


10 
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30 
0.087 
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40 
10 
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1 


5 
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0.30 

9 X 10"^ 
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0.05 


65 
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TABLE I: Nominal parameters for single-photon generation with bare diamond optical centers [indicated by o] and cavity- 
enhanced centers [(•)] with V = (Ac/n)^. Pulsed-excitation parameters are chosen so that multi-photon emission probability 
Pm — O(10~^). * denotes the estimated dipole moment for the zero-phonon transition and f the total dipole moment of the 
center across all phonon-assisted transitions. 



and implementation requirements of these centers, in- 
cluding the NV center from Ref. in Table H The 
assumed setup of (T, r) = (2.3 ps, 1.3 THz) implies a mod- 
est peak intensity 1 GW/cm^ and pulse energy of 25 pJ 
on the NV, Pi = 0.95, and P^ = 0(10-^). Pulses of 6 pJ 
and T = 0.6 ps would give Pi = 0.54 and P^ = 0(10"'^). 
We have also included the performance of the bare cen- 
ters in the tabulation and their emission probabilities are 
estimated using Pi = 27r/(r7totai)- 



III. QUANTUM KEY DISTRIBUTION 

Robust and high-performance SPSs are an invaluable 
resource for implementing photon-based QKD protocols 
such as the BB84 [8]. In this section, we are interested 
in exploring the potential of diamond-based sources in 
a range of implementations by using BB84 as a met- 
ric for comparisons. Minimizing multi-photon emission 
probability P^ is fundamentally important for extend- 
ing the range of secure communication, we therefore as- 
sume that the cavity-enhanced diamond centers are ex- 
cited with suitable excitation (T, r) (see Sec. II) whereby 
Pm = O(10~^) at the cost of reduced Pi ~ 0.6, and 
without further attenuation, g^'^^ = O(10~^). For the 
bare centers, we let g^'^^ = 0. 

QKD allows two parties, Alice (sender) and Bob (re- 
ceiver) to establish a secret key via a quantum channel 
that is open to an eavesdropper. Alice uses single-photon 
states to communicate classical information by encoding 
each photon in one of four polarizations: 0°, 45°, 90°, 
135° orientations. The first two correspond to classical 
bit while the latter two bit 1. Alice (Bob) randomly 
chooses one of two bases, rectilinear or diagonal, to pre- 
pare (measure) the photons. Consequently, they commu- 
nicate openly the bases they have used and retain those 
bits with matching bases. Finally, both parties perform 
local classical post-processing, namely error correction 
and privacy amplification, to correct erroneous bits and 
reduce the bit-string to a shared secure key asymptoti- 



cally unknown to the eavesdropper. 

In the fiber-based scheme, time-bin encoding of signal 
states is preferred because single-mode fiber is susceptible 
to polarization loss. In contrast, polarization encoding is 
suitable for free-space QKD owing to the non-birefringent 
nature of the atmosphere. In both cases, the secure key 
generation rate Q (bits s~^) of the protocol using a sub- 
Poissonian source is treated by Waks et al. in Ref. [6^, 

^sps = vqPi [/3T(e) - f{e)h{e)] (7) 

where v is the repetition rate for bitstream single photons 
and e is the total bit error rate, q is the fraction of the 
correctly measured signals (typically 1/2 because half the 
time Alice and Bob disagree on the bases chosen). The 
compression function r(e) of the privacy amplification 
procedure is 

T(e) = -log2 [1/2 + 2e//3 - 2{e/pf] . (8) 

h{e) is the Shannon entropy function of a bit. Pd is 
the probability that Bob registers a detection event for 
a signal pulse. Although the mean number n of photons 
emitted from Alice's source is (Pi + 2Pm) ^ Pi, she may 
add extra attenuation ^ to reduce the relative number 
of multi-photon states. As a result, the key rate can be 
optimized with respect to this adjustable parameter ^ for 
a given communication range. Therefore, the detection 
probability at Bob is Pd = (,r]Pi-\-N, where r] is the total 
optical loss in the system, N is the total noise from the 
channel, the environment and the detection units. 

In Eq. [SI is the fraction of detection events origi- 
nating from single-photon states. Expressing in terms 
of the attenuation- free parameter the multi-photon 
emission probability becomes Pm = ^^Pf^^^V^ and (3 = 
(Pd — Pm)/Pd [68]. The function /(e) characterizes the 
performance of the error correction algorithm for a bit 
error rate e = (eo^/^Pi + A^/2)/Pd, where the baseline 
error rate eo (typically 2%) accounts for erroneous de- 
tection and state preparation. Signals can be lost dur- 
ing channel and apparatus transmission or coupling into 



6 



and out of the channel, and due to detection inefficiency: 
V = VoVcViVd^ where rjo (typically 60%) is the total opti- 
cal transmittance at Bob's end, r]c the total coupling effi- 
ciency, r]i the channel transmittance, and rjd the quantum 
efficiency of a detector. 

For a quasi-SPS, Gsps scales with loss as r]. When the 
upper bound on the allowable error rate is taken to be 
(3/4:, the loss cutoff for a sub-Poissonian source is [68] 

1 /TV Pig^^h 

In other words, key generation between the parties re- 
mains possible if 7^ > r^min so that N and g'^'^^ should be 
minimized for long-range communication. 

On the other hand, the corresponding key rate formula 
for a weak coherent source (WCS) is [69| 

Gwcs = i^qPd [pr{e/(3) - f{e)h{e)] (10) 

where for a coherent state with mean photon number n, 
Pd = 1 — exp(— 7772) -\-N when r] is small and noise level is 
the same for all photon states. The error rate e = {eo[l — 
ex.-p{r]n)] + N/2}/P^ and the probability of emitting a 
multi-photon state from Alice is Pm = 1 — (1+n) exp(— n). 
Since the secure key is still extracted from the single- 
photon states, the source must be appropriately attenu- 
ated prior to transmission and the optimal photon num- 
ber that maximizes the key rate is nopt ^ ^, leading to 
the scaling 5wcs = 0{rf). The estimate for the corre- 
sponding loss cutoff is T^min = V^27V /(I — 4eo) [70]. 

In the decoy state QKD, Alice uses additional vacua 
and multi-photon states to estimate the channel loss and 
error rate for each n-photon signal accurately [TO, Ei so 
that extra attenuation of the source is not required. The 
optimal photon number nop 0.7 for eo = 2% so that 
the key rate is now of 0{r]) 

awcs,decoy = ^^Pd {/?[! - Ke/[3)] - f{e)h{e)} . (11) 

The noise level and optical loss vary greatly with spe- 
cific QKD realizations. Here we explore three scenarios, 
namely QKD over (A) a silica fiber-optic telecom chan- 
nel, (B) a free-space terrestrial channel, and (C) satellite- 
ground uplink and downlink. In each case we compare 
the performance realized using cavity-enhanced and bare 
diamond centers, with WCS counterparts at wavelengths 
650 nm and 1.55 /im. 

The decoy state QKD protocol through a fiber link 
has been demonstrated by several groups. Rosenberg et 
al. achieved 107 km range at 1.3 /im wavelength using 
high-performance transition-edge sensor (TES) photode- 
tectors [7l|. Using the parameters in Ref. exchange 
over 140 km at 1.55 /im is feasible [ll^ with dark count 
rate of 100 Hz and 3.5 ns gate time. Dixon et al. achieved 
1 Mbit s~^ (10.1 kbit s~^) secure key rate over a telecom- 
fiber distance of 20 km (100 km) using an InGaAs photon 
detector at the same wavelength [73]. Single NV centers 
in diamond nanocrystals have been used to realize se- 
cure exchange over short distances of 30-50 m in free 



space [iH, |26|. Long-distance exchange (over 144 km 
at 850 nm) has only been achieved using attenuated 
pulses and decoy states [74]. However, as we will show, 
while the decoy state protocol with high-speed WCSs are 
more suitable for practical QKD implementations than 
bare diamond centers, cavity-enhanced centers can de- 
liver comparable performance and in some cases, outper- 
form WCS-based implementations. We note that the fea- 
sibility and engineering of ground-satellite QKD system 
have been investigated by several authors [75|, [t^, [73, [7^1 • 
The first experimental study of such a setup has also been 
reported in Ref. [79|. 



A. Silica-fiber based QKD 

In the fiber-based QKD scheme, the noise is largely 
due to dark counts at Bob's detectors. In a typical four- 
detector implementation, N = ARAt where R is the dark 
count rate and At is the size of measurement gate win- 
dow, usually chosen to match photon pulse length 
Low dark counts have been achieved in various detec- 
tion implementations. Widely available silicon avalanche 
photodiodes have a dark count rate P ~ 25 counts s~^ 
and ~ 65% efficiency over the visible regime [80]. De- 
tectors well-suited for near-IR wavelengths are supercon- 
ducting nanowire single-photon detectors (SSPD) [8l[ 
and TES [82|. At 1.55 /im, a TES has been shown to 
deliver 50% quantum efficiency for an erroneous back- 
ground count rate of 0.053 Hz whereas an SSPD can de- 
liver 57% efficiency with higher dark counts ~ 1000 Hz. 

For practical distances, the dominant loss is due to 
fiber attenuation, given by r]i = io~"^/^^ where a 
is characteristic attenuation coefficient and / is fiber 
length. The choice of silica-based optics favors the near- 
IR wavelength with a = 0.2 dB km~^, compared to 
2.5 — 6 dB km~^ in the visible. The bare centers suf- 
fer from low coupling efficiency into fiber modes such 
that r]c ^ 1% for single photons. The directionality and 
coherence of the photons emitted via cavity mode can 
improve this efficiency to over 50% [83]. For the WCS, 
our choice of n = nopt implies that the combined atten- 
uation at Alice's end has already included WCS-to-fiber 
coupling loss so that r]c for WCS can be set to unity. 

We summarize the key parameters and performance for 
the fiber-based systems with different sources in Table HIl 
Equations [THTTl are used to calculate the key rates as a 
function of total optical loss r^dB = — lOlogio^ (i-^- 
decibel) in Fig.[4l The curves corresponding to WCSs are 
plotted with their respective optimal photon numbers. 
Here we compare transmissions clocked at the maximum 
repetition rate u of the diamond-based SPSs with WCSs 
driven at rate 10 GHz. In practice, the repetition rate 
is limited by the timing jitter (< 0.1 ns) and reset time 
(< 1 /is) of the detectors but this effect can be ignored 
when these time scales are shorter than the mean arrival 
times of photons after a lossy channel. The prescrip- 
tions used here are still valid if any signals arrived during 
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Source 




At 

(ns) 


a 

(dB/km) 




N 

(1 /pulse) 


'^min 

(dB) 


NV 


(•) 

o 


0.5 

420 


6 


0.5 
0.01 


5.0 X 10"^ 
4.2 X 10"^ 


68 
28 


NE8 


(•) 

o 


0.03 
16.4 


2.5 


0.5 
0.01 


3.3 X 10-^ 
1.6 X 10"^ 


74 

55 


SiV 


(•) 

o 


0.02 
3.4 


3.5 


0.5 
0.01 


1.5 X 10"^ 
3.4 X 10"'^ 


75 
50 


WCS* 
WCS^ 




0.1 


6 

0.2 


1 
1 


1.0 X 10-^ 

2.1 X 10"^^ 


66 
92 



TABLE II: Nominal parameters for BB84, based on silica- 
fiber that is optimized for 1.55 /x, with different single-photon 
sources. The bare centers are indicated by o and cavity- 
enhanced centers (•). * denotes 10 GHz WCS at 650 nm 
and t at 1.55 /xm. The values of loss cutoff ?7min for WCSs 
correspond to decoy state implementations. 



Silica fiber-based QKD 

— Cavity-SiV 
-- SiV 

— Cavity-NE8 
--NE8 

— Cavity-NV 
--NV 

•••• WCS 1.55 ^im 
■ - (decoy state) 
•••• WCS 650 nm 
.. , , • - (decoy state) 
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FIG. 4: (Color online) Secure key rate over silica- fiber link 
as a function of total loss, for different single-photon sources. 
Dashed: BB84 with bare optical centers NV (black), NE8 
(blue), and SiV (red). Solid: BB84 with cavity-enhanced cen- 
ters distinguished by same color curves. Dotted: BB84 with 
10 GHz WCSs at 650 nm (green) and 1.55 /xm (pink). Dashed- 
dotted: Decoy state with WCSs distinguished by same color 
curves. We use the parameters from Table Ull rja — 0.65 (0.5), 
R — 25 Hz (0.035 Hz) for visible (near-IR) wavelength trans- 
mission, rjo = 0.6 and eo = 0.02 for all cases. To aid compar- 
isons, solid circles are used to indicate the key rates over fiber 
length / of 5 km. 



the reset time are discarded. To avoid fiber dispersion, 
only the emissions via ZPL transition of the bare centers 
are considered for QKD, and the ZPL spontaneous decay 
rates of the bare centers are estimated to be 2.5 x 27r MHz 
(NV), 60 X 27r MHz (NE8) and 300 x 27r MHz (SiV). 

We obtain substantial improvements in terms of key 
rate (by 3-4 orders of magnitude) and range by using 
cavities. Over a 5 km fiber, the NE8 center delivers 
62 Mbits s-^ key rate, the SiV 7 Mbits s"^ and the 
NV 0.08 Mbits s~^. However, since the telecom fiber 
favors near-IR transmission, a 1.55 jam WCS with decoy 
states is expected to be superior with over 200 Mbits s~^ 
key rate for / = 5 km. Given an ultralow noise level of 



(a) m 



^r^^ horizontal path / -^^'^ 
^^ Jr turbulent atm ^ , 0-1^^-2/3 4 
alt. 3km *^«~^^ ^ ^ 



(c) Alice 



LEQ,. « 
satellite 4 

(alt! H~2000km) 




FIG. 5: (Color online) Generic (a) open-air terrestrial and 
(b) LEO satellite-ground communications over a path length 
/ through atmospheric turbulence characterized by refractive- 
index structure parameter C^- (c) Transmitter using a typical 
Cassegrain telescope. Single photons generated at a repetition 
rate u are encoded using electro-optic modulator, (d) Receiv- 
ing telescope with an angular field of view ^fov sees a broad- 
ened spot with radius WeH- Measurement basis is randomly 
chosen with a beamsplitter (BS) and polarization states are 
distinguished with polarizing beamsplitters (PBSs). ATP sys- 
tem, post-processing, key management are not depicted here. 



2.1 X 10~^^ pulse"^ (i.e . 5 orders of magnitude less than 
that used in Refs. [h], ^), we estimate an optimistic 
420 km communication range using the decoy states. 



B. Open-air terrestrial QKD 

An open-air system is susceptible to stray photons 
and, in the visible and near-IR regimes, sunlight is 
the dominant source of such stray photons during the 
day, and starlight or artificial light sources dominate at 
night. The background brightness i^b observed from the 
ground varies from 1.5 Wm~^Sr~^/im~^ on a clear day- 
time to 1.5 X 10~^Wm~^Sr~^/im~^ on a full moon clear 
night [t^I- However, the use of time-gate, frequency and 
spatial filters, together with an active high-precision ac- 
quisition, tracking and pointing (ATP) system between 
the transmitting (Alice) and receiving (Bob) telescopes, 
should improve the signal-to-noise ratio and allow com- 
munication over long distances [77]. Solar background 
can also be reduced by operating in the Fraunhofer win- 
dows. A setup is schematically depicted in Fig. [5l The 
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gated noise counts per pulse is then given by [77 

N = (iJbf^fov^rec-Bfilter + 4i?)At, 



(12) 



where Bob uses a receiving telescope with a field of view 
(FOV) r^fov, a numerical aperture of diameter Dn and 
area /irec = 7TDj^/4, gate sizes (temporal filter) of At 
and frequency filters of width ^filter- For a detector of 
diameter and a telescope with a focal length / to 
aperture size ratio of a = f /Dr, the angular FOV 6fov = 
Dd/{DRa) and the solid angle l^fov = 27r(l — cos^fov) ~ 
tt^Iq^ for a large / ^ D^. When the telescope is used 
with an eyepiece, 6>fov is the ratio of the diameter of the 
field stop to the focal length, as illustrated in Fig. [SJd). 

Signal attenuation through the atmosphere is at- 
tributed to scattering by aerosols and absorption by at- 
mospheric gases. In the visible and near-IR, the atmo- 
sphere is generally transmissive and the major gaseous 
absorber is water vapor. In particular, there are low- 
loss windows at around 770 and 870 nm The typical 
values for attenuation by scattering at sea level can be 
found in Ref. [84] for a range of weather conditions and 
visibility, e.g., a - 0.4 - 0.6 dB km"^ (- 2.2 dB km"^) 
for a visibility range of 23.5 km (5 km) in a standard- 
clear (medium-haze) condition. These values decrease 
with altitude, e.g., by one order of magnitude at ~ 3 km 
above sea level. To estimate spectral transmittance of 
H2O, one can refer to tabulations of experimentally de- 
termined transmission data (e.g. Ref. or the at- 
mospheric model in the MODTRAN code [86|]. As our 
reference we take the values corresponding to tempera- 
ture 25° C and relative humidity 60%, the absorption loss 
is roughly —12 dB in the visible and —4 dB at 1.55 /im 
over a 150 km path at 3 km altitude. 

The signal beam leaving the transmitting telescope 
diffracts and diverges to give rise to pulse spreading. For 
a Gaussian beam of wavelength A, the divergence half- 
angle is = X/ (ttil'o), where the radius of the beam waist 
Wo should be narrower than half of the diameter Dt of 
the collimating telescope (we let wq = 0.35Dt) so that 
most of the beam is transmitted. When the signal pho- 
tons are detected at a distance / away. Bob sees a broad- 
ened spot with diffraction-limited radius (~ 01). Fur- 
thermore, atmospheric turbulence induces random dis- 
placement of the beam centroid (beam wander) when 
the beam propagates through optical turbules that are 
smaller than the beam diameter. This effect leads to a 
long-term beam radius [STj 



(13) 



where k = 27r/A is the wave number and p = 
[0.55/c^/C^]~^/^ for a fiat profile of the refractive index 
structure along the propagation path. The typical val- 
ues of refractive-index structure parameter C^, which 
measures the strength of fiuctuation in n, range from 
-13 ^j^der weak to stronger turbulent 



10" 



-17 
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conditions. Thus the fraction of the light collected at the 
receiving telescope is 



(b) 



(c) 



Source 


N ?7min 

(1/pulse) (dB) 


N ?7min 

(1/pulse) (dB) 


N ?7min 

(1/pulse) (dB) 


NV 




6.0 X IQ-s 68 
5.0 X 10"^ 27 


5.6 X 10-^ 68 

4.7 X 10"^ 28 


6.1 X 10-^ 58 
5.1 X 10"^ 17 


NE8 




1.6 X IQ-^ 72 
2.0 X 10"^ 54 


1.1 X 10"^ 73 
1.9 X 10-^ 55 


7.4 X 10"^ 58 

2.5 X 10"^ 43 


siv 




1.2 X 10-^ 73 
4.2 X 10"^ 49 


7.9 X 10-^ 74 
3.9 X 10-^ 50 


6.3 X 10-^ 58 
4.8 X 10"^ 38 


wcs* 
wcs^ 


1.2 X 10"^ 65 
4.8 X 10"^ 69 


1.1 X 10"^ 66 
2.9 X 10"^ 71 


1.2 X 10"^ 54 
2.9 X 10"'^ 51 



Vc 



1 



-^H/(2<f) 



(14) 



TABLE III: Estimated noise N and loss cutoff ?7min for (a) 
open-air terrestrial at 3 km altitude, (b) LEO (low earth or- 
bit) satellite-ground uplink and (c) downlink implementation 
of the BB84 QKD protocol with different photon sources. We 
assume a night-time operation where the background bright- 
ness ifb = 1-5 X 10~^Wm~^Sr~Vrn~^- The values of ?7min for 
WCSs correspond to decoy state implementations. 



In the open-air terrestrial QKD setup, we consider a 
link at ~ 3 km above sea- level on a standard-clear, full 
moon night (i^b = 1-5 x 10~^ Wm~^Sr~^/im~^, visibil- 
ity ~24 km). Following the Canary islands experiment 
reported in Ref. [t^ , the emitted photons are collimated 
by 0.15 m optics and then intercepted by a 1 m diame- 
ter, //39 telescope (39 m focal length) that has a FOV 
of 12.5 /irad for a typical detector size 0.5 mm. We es- 
timate a moderate ~ 4 x 10~^^ m~^/^ that causes a 
roughly —14 dB collection loss (i.e. Wef[ ~ 3.5 m) in that 
demonstration. Under these conditions, the magnitude 
of the beam wander is of same order as the diffraction- 
limited spot size at about 10 km (40 km) for the visible 
(1.55 /im) wavelengths, and beyond which beam- wander 
effect to dominate. Thus, while shorter wavelengths suf- 
fer less diffraction, this advantage diminishes for long 
distances > 60 km unless Cn improves. For a typical pi- 
cosecond time-of- arrival jitter, the setup can be clocked 
at the proposed repetition rates. Given the availabil- 
ity of high-transmittance (> 90%) atom filters of width 
0.01 nm, we set ^filter = 0.01 nm or AA of sources with 
wider spectral widths. 

The estimated loss cutoffs of this setup with differ- 
ent sources are tabulated in Table llllt and the key 
rates are plotted in Fig. [D^a). Over a 150 km link, 
NE8 and SiV centers can deliver a remarkable 0.8 and 
1.8 Mbits s~^ respectively, whereas 10 GHz 650 nm 
and 1.55 jam WCSs with decoy states enables 0.1 and 
1.2 Mbits s~^. The total losses of around —35 dB in- 
clude an average r]c = —14.5 dB due to beam spread- 
ing, — 7 dB due to atmospheric loss via scattering, —4 
(near-IR) to —12 dB (visible) from H2O absorption, and 
—4.5 dB due to apparatus imperfections. Since the loss 
cutoff ?7min is arouud —70 dB during the night time, we 
predict over 300 km exchange is possible. On a clear 
day (i^b = 1.5 Wm-^Sr" Vm"^), r^min ^ -46 dB would 



9 




^dB [dB] 77dB [dB] T]dB [dB] 



FIG. 6: (Color online) QKD performance of (a) open-air, terrestrial link at 3 km above sea-level, (b) free space uplink between a 
transmitter at an altitude ho and a LEO satellite at an altitude 2000 km, and (c) LEO satellite-ground downlink, on a standard- 
clear night with brightness Hh — 1.5 x 10~^Wm~^SY~^ /jm~^ . In (a), Dt — 0.15 m at Alice, Dr = 1 m, Dd — 0.5 mm, a — 39, 
and FOV = 12.5 /xrad at Bob. In (b), = 1 m on the ground (Alice), Dr = 0.1 m, Dd = 0.5 mm, a = 50 and FOV = 
100 /xrad on the satellite (Bob). In (c), Dt = 0.1 m on the satellite, Dr = 1 m, Dd = 0.5 mm, a = 5 and FOV = 100 /irad on 
the ground. The values for r/d, R, Vo and eo follow Fig. 2] noise rates N are tabulated in Table Hill Circles are used to indicate 
the estimated key rates for a particular transmission distance or altitude. 



still allow exchanges over distances < 150 km with the 
diamond-based sources. 



C. Satellite-ground QKD 

Under a standard clear-sky condition, transmission 
through the entire atmosphere between a ground-based 
station at an altitude ho and a satellite H = h^^l can 
be as high as 60-70% at the visible and 85% for 1.55 /im. 
Signals are largely lost at collection due to pulse spread- 
ing and beam wander. Along a zero-zenith angle path, 
the effective spot size is [88] 

^eff = ei^i^7.7b^{k/e^Y/\ (15) 

where ji is related to varying strength of optical turbu- 
lence C^{h) as a function of altitude h. For uplink paths, 

while for downlink paths, 

and according to the Hufnagel- Valley atmospheric model, 

Clih) = 0.00594(w/27)2(10-^/i)^°e-''/i°™ + 

2.7 X 10-1 V^/i^™ + C2(0)e-'^/i™ (18) 

where h is in meters, v is the pseudowind speed that 
ranges from 10 - 30 m s~^ for moderate to strong wind, 
and nominal value of ground turbulence is C^(0) = 
1.7 X 10~^^ m~^/^. For satellite uplinks, beam wander is 
significant as the beam diameter is smaller than many of 



the turbules encountered. On the other hand, downlinks 
are less affected as the signal beam is already broadened 
by diffraction before entering the Earth's atmosphere. 

Considering a satellite's allowed payload, we assume 
the use of 10 cm optics on the satellite [76]. In the 
uplink communication between a low earth orbit (LEO) 
satellite at = 2000 km and a ground-based telescope 
ho = [Fig. EJb)], Alice uses aim diameter telescope 
{wo = 0.35 m) to transmit signal photons to //50 down- 
ward pointing optics on the satellite with a 0.5 mm de- 
tector. This setup implies an effective beam spot radius 
of roughly 31 m for both visible and near-IR wavelengths, 
and collection loss (Eq.[H|) of —53 dB. In general, Wef[ de- 
creases rapidly with increasing transmitter beam radius 
for small beams {wo < 4 cm) but this reduction becomes 
less pronounced for larger apertures {wo > 4 cm). The 
collection loss can be reduced significantly by situating 
Alice at higher altitudes, e.g., ho = 3 km yields ~3 m 
beam at = 2000 km and ~ —30 dB loss. Since the 
receiving optics has a wide FOV of ^fov = 100 /irad, re- 
quirements on the tracking accuracy between the satellite 
and the ground station are not stringent. 

Figure[6](b) shows the key rates for operation on a stan- 
dard clear night. While the setup is not optimized, we 
see that the centers can deliver up to 24 kbits s~^ key 
rate from the ground and over 3 Mbits s~^ from 3 km 
altitude, which are one magnitude faster than 10 GHz 
WCSs. The cavity-enhanced sources and higher altitude 
implementations also allow operations even in hazy day- 
time condition with 15 Wm~^Sr~^/im~^ brightness and 
^min ^ -40 dB. 

For the downlink case, the signal beam leaves the satel- 
lite with a spot radius of 3.5 cm and is taken to be 
detected by a //5 1 m telescope, where Wef^ = 12 (for 
637 nm) to 28 m (1.55 /im) spot radius, or —24 and 
—32 dB losses respectively. Since the turburlence is much 
weaker at high altitudes, we expect from Eqs. [T5land[T71 
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that the spot size is essentiahy the same as the diflFrac- 
tive spot size iL'eff ^ 01 at arrival. Thus the cohection loss 
only improves significantly by using larger transmitting 
and receiving apertures, as opposed to using higher al- 
titude receivers. Notably, it is advantageous to perform 
downlink at visible wavelengths because the divergence 
angle is larger at 1.55 /im, and this leads to a lower collec- 
tion loss by near 8 dB. Fig. ^c) shows that the centers 
can produce up to 0.6 — 13 Mbits s~^ key rates. Sim- 
ilarly, downlink can also be performed on a clear day 
with 77min ~ — 29 dB for the enhanced centers. Finally 
we summarize the expected noise levels and loss cutoffs 
for satellite-ground QKD at night time in Table Hill 

IV. CONCLUSION 

We have studied the effect of a cavity on diamond 
defects, namely an NE8 and SiV optical centers, in en- 
hancing its spectral properties for high-performance sin- 
gle photon generation. The Purcell enhancement in the 
high-Q limit [Q = O(IO^), V = 0{X^)] of the weak cou- 
pling regime enables rapid photon generation on picosec- 
ond time scale. By applying suitable excitation schemes, 
these systems constitute high performance and fidelity 
single-photon sources that are capable of operating at 
> 10 GHz repetition rates. 

Standard BB84 using true single-photon sources pro- 
vides an important benchmark implementation for quan- 



tum key distribution. In this work, we have contrasted 
bare and cavity-enhanced centers (including the NV cen- 
ter) with the weak coherent states in various prospec- 
tive long-distance fiber-based and open-air setups. Fiber- 
optical communications with these diamond centers are 
expected to be short-ranged due to high loss through IR- 
optimized silica fiber at visible wavelengths. However, 
in open-air arrangements, cavity-centers can deliver per- 
formance that is comparable or better than the decoy 
state implementation. Notably, during day and night- 
time operation, these centers can achieve over 1 Mbits s~^ 
secure key rate in low earth orbit satellite-ground ex- 
changes, and terrestrial point-to-point communications 
over 300 km path at 3 km altitude. 

Acknowledgments 

We thank W. J. Munro, C. T. Chantler, A. Roberts, Z. 
W. E. Evans, A. M. Stephens, D. A. Simpson and B. C. 
Gibson for helpful discussions, and F. M. Hossain for the 
figures of the NE8 and SiV structures in Figs.[T](c,d). We 
acknowledge the support of Quantum Communications 
Victoria, funded by the Victorian Science, Technology 
and Innovation (STI) initiative, the Australian Research 
Council (ARC), and the International Science Linkages 
program. A.D.G. and L.C.L.H. acknowledge the ARC 
for financial support (Projects No. DP0880466 and No. 
DP0770715, respectively). 



[1] D. F. Walls and G. J. Milburn, Quantum Optics 

(Springer-Verlag, New York, 1995). 
[2] E. S. Polzik, J. Carri, and H. J. Kimble, Phys. Rev. Lett. 

68, 3020 (1992). 
[3] V. Giovannetti, S. Lloyd, and L. Maccone, Science 306, 

1330 (2004). 

[4] P. G. Kwiat, K. Mattle, H. Weinfurter, A. Zeilinger, A. 
V. Sergienko, and Y. H. Shih, Phys. Rev. Lett. 75, 4337 
(1995). 

[5] L.-M. Duan, M. Lukin, J. Cirac, and P. Zoller, Nature 
(London) 414, 413 (2001). 

[6] E. Knill, R. Laflamme, and G. J. Milburn, Nature (Lon- 
don) 409, 46 (2001). 

[7] N. Gisin, G. Ribordy, W. Tittel, and H. Zbinden, Rev. 
Mod. Phys. 74, 145 (2002). 

[8] C. H. Bennett and G. Brassard, in Proceedings of IEEE 
International Conference on Computers, Systems and 
Signal Processing (IEEE, New York, 1984), p. 175. 

[9] W.-Y. Hwang, Phys. Rev. Lett. 91, 057901 (2003). 
[10] X.-B. Wang, Phys. Rev. Lett. 94, 230503 (2005). 
[11] H.-K. Lo, X. Ma, and K. Chen, Phys. Rev. Lett. 94, 
230504 (2005). 

[12] T. Horikiri and T. Kobayashi, Phys. Rev. A 73, 032331 
(2006). 

[13] C. K. Hong and L. Mandel, Phys. Rev. Lett. 56, 58 
(1986). 

[14] J. Kim, O. Benson, H. Kan, and Y. Yamamoto, Nature 
(London) 397, 500 (1999). 



[15] B. T. H. Varcoe, S. Brattke, M. Weidinger, and H. 

Walther, Nature (London) 403, 743 (2000). 
[16] M. Hijlkema, B. Weber, H. P. Specht, S. C. Webster, A. 

Kuhn, and G. Rempe, Nat. Phys. 3, 253 (2007). 
[17] M. Keller, B. Lange, K. Hayasaka, W. Lange, and H. 

Walther, Nature (London) 431, 1075 (2004). 
[18] C. Brunei, B. Lounis, Ph. Tamarat, and M. Orrit, Phys. 

Rev. Lett. 83, 2722 (1999). 
[19] B. Lounis and W. E. Moerner, Nature (London) 407, 

491 (2000). 

[20] A. J. Shields, Nat. Photonics 1, 215 (2007). 

[21] C. Kurtsiefer, S. Mayer, P. Zarda, and H. Weinfurter, 

Phys. Rev. Lett. 85, 290 (2000). 
[22] T. Gaebel, L Popa, A. Gruber, M. Domhan, F. Jelezko, 

and J. Wrachtrup, New J. Phys. 6, 98 (2004). 
[23] E. Wu, J. R. Rabeau, G. Roger, F. Treussart, H. Zeng, 

P. Grangier, S. Prawer, and J.-F. Roch, New J. Phys. 9, 

434 (2007). 

[24] C. Wang, C. Kurtsiefer, H. Weinfurter, and B. Burchard, 
J. Phys. B: At. Mol. Opt. Phys. 39, 37 (2006). 

[25] A. Beveratos, R. Brouri, T. Gacoin, A. Villing, J.-P. 
Poizat, and P. Grangier, Phys. Rev. Lett. 89, 187901 
(2002). 

[26] A. Alleaume, F. Treussart, G. Messin, Y. Dumeige, J.-F. 

Roch, A. Beveratos, R. Brouri- Tualle, J.-P. Poizat, and 

P. Grangier, New J. Phys. 6, 92 (2004). 
[27] S. Kiihn, C. Hettich, C. Schmitt, J.-PH. Poizat, and V. 

Sandoghdar, J. Microsc. 202, 2 (2001). 



11 



[28] J. R. Rabeau, S. T. Huntington, A. D. Greentree, and S. 

Prawer, AppL Phys. Lett. 86, 134104 (2005). 
[29] S. Trpkovski, L. C. L. Hollenberg, and S. Prawer, Laser 

Focus World 44, 70 (2008). 
[30] P. Neumann, R. Kolesov, V. Jacques, J. Beck, J. Tisler, 

A. Batalov, L. Rogers, N. B. Manson, G. Balasubrama- 

nian, F. Jelezko, and J. Wrachtrup, New J. Phys. 11, 

013017 (2009). 

[31] T. Feng and B. D. Schwartz, J. Appl. Phys. 73, 1415 
(1993). 

[32] L Aharonovich, C. Zhou, A. Stacey, J. Orwa, S. Castel- 
letto, D. Simpson, A. D. Greentree, F. Treussart, J.-F. 
Roch, and S. Prawer, Phys. Rev. B 79, 235316 (2009). 

[33] D. A. Simpson, E. Ampem-Lassen, B. C. Gibson, S. Trp- 
kovski, F. M. Hossain, S. T. Huntington, A. D. Greentree, 
L. C. L. Hollenberg, and S. Prawer, Appl. Phys. Lett. 94, 
203107 (2009). 

[34] L Aharonovich, S. Castelletto, D. A. Simpson, A. Stacey, 
J. McCallum, A. D. Greentree, and S. Prawer, Nano Lett. 
9, 3191 (2009). 

[35] C.-H. Su, A. D. Greentree, and L. C. L. Hollenberg, Opt. 

Express 16, 6240 (2008). 
[36] S. Noda, M. Fujita, and T. Asano, Nat. Photonics 1, 449 

(2007). 

[37] S. Tomljenovic-Hanic, M. J. Steel, C. Martijn de Sterke, 
and J. Salzman, Opt. Express 14, 3556 (2006). 

[38] L Bayn and J. Salzman, Eur. Phys. J. Appl. Phys. 37, 
19 (2007). 

[39] C. Kreuzer, J. Riedrich-Moller, E. Neu, and C. Becher, 

Opt. Express 16, 1632 (2008). 
[40] L Bayn and J. Salzman, Opt. Express 16, 4972 (2008). 
[41] S. Tomljenovic-Hanic, A. D. Greentree, C. Martijn de 

Sterke, and S. Prawer, Opt. Express 17, 6465 (2009). 
[42] P. E. Barclay, K.-M. Fu, C. Santori, and R. G. Beausoleil, 

Opt. Expr. 17, 9588 (2009). 
[43] A. D. Greentree, B. A. Fair child, F. Hossain, and S. 

Prawer, Mater. Today 11, 22 (2008). 
[44] C. F. Wang, R. Hanson, D. D. Awschalom, E. L. Hu, T. 

Feygelson, J. Yang, and J. E. Butler, Appl. Phys. Lett. 

91, 201112 (2007). 
[45] M. Barth, N. Niisse, B. Lochel, and O. Benson, Opt. Lett. 

34, 1108 (2009). 
[46] P. E. Barclay, K.-M. C. Fu, C. Santori, and R. G. Beau- 
soleil, e-print arXiv:0908.2148 (2009). 
[47] V. R. Almeida, Q. Xu, C. A. Barrios, and M. Lipson, 

Opt. Lett. 29,1209 (2004). 
[48] J. T. Robinson, C. Manolatou, L. Chen, and M. Lipson, 

Phys. Rev. Lett. 95, 143901 (2005). 
[49] M. P. Hiscocks, C.-H. Su, B. C. Gibson, A. D. Greentree, 

L. C. L. Hollenberg, and F. Ladouceur, Opt. Express 17, 

7295 (2009). 

[50] C. D. Clark, H. Kanda, L Kiflawi, and G. Sittas, Phys. 

Rev. B 51, 16681 (1995). 
[51] A. V. Turukhin, C. H. Liu, A. A. Gorokhovsky, R. R. 

Alfano, and W. Phillips, Phys. Rev. B 54, pp. 16448 

(1996) . 

[52] C. K. Law and H. J. Kimble, J. Mod. Opt. 44, 2067 

(1997) . 

[53] E. M. Purcell, Phys. Rev. 69, 681 (1946). 

[54] S. Haroche and J. M. Raimond, in Advances in Atomic 

and Molecular Physics (Academic Press, New York, 

1985), Vol. 20, p. 350. 
[55] G. Cui and M. Raymer, Opt. Express 13, 9660 (2005). 
[56] Y. Dumeige, F. Treussart, R. Alleaume, T. Gacoin, J.-F. 



Roch, P. Grangier, J. Lumin. 109, 61 (2004). 
[57] T. Brabec and F. Krausz, Rev. Mod. Phys. 72, 545 
(2000). 

[58] A. Baltuska and T. Kobayashi, Appl. Phys. B 75, 427 

(2002) . 

[59] U. Keller, Nature (London) 424, 831 (2003). 

[60] G. Cerullo and S. De Silvestri, Rev. Sci. Instrum. 74, 1 

(2003) . 

[61] A. Bartels, D. Heinecke, and S. A. Diddams, Opt. Lett. 

33, 1905 (2008). 
[62] A. D. Greentree, J. Salzman, S. Prawer, and L. C. L. 

Hollenberg, Phys. Rev. A 73, 013818 (2006). 
[63] M. J. Fernee, H. Rubinsztein-Dunlop, and G. J. Milburn, 

Phys. Rev. A 75, 043815 (2007). 
[64] C.-H. Su, A. D. Greentree, W. J. Munro, K. Nemoto, and 

L. C. L. Hollenberg, Phys. Rev. A 78, 062336 (2008). 
[65] C.-H. Su, A. D. Greentree, W. J. Munro, K. Nemoto, and 

L. C. L. Hollenberg, Phys. Rev. A 80, 033811 (2009). 
[66] L. Tian and H. J. Carmichael, Phys. Rev. A 46, R6801 

(1992). 

[67] H. J. Carmichael, An Open System Approach to Quantum 

Optics (Springer, Berlin, 1993). 
[68] E. Waks, C. Santori, and Y. Yamamoto, Phys. Rev. A 

66, 042315 (2002). 
[69] D. Gottesman, H.-K. Lo, N. Liitkenhaus, and J. Preskill, 

Quant. Inf. Comp. 5, 325 (2004). 
[70] G. Brassard, N. Liitkenhaus, T. Mor, and B. C. Sanders 

Phys. Rev. Lett. 85, 1330 (2000). 
[71] D. Rosenberg, J. W. Harrington, P. R. Rice, P. A. 

Hiskett, C. G. Peterson, R. J. Hughes, A. E. Lita, S. W. 

Nam, and J. E. Nordholt, Phys. Rev. Lett. 98, 010503 

(2007) . 

[72] C. Cobby, Z. L. Yuan, and A. J. Shields, Appl. Phys. 

Lett. 84, 3762 (2004). 
[73] A. R. Dixon, Z. L. Yuan, J. F. Dynes, A. W. Sharpe, and 

A. J. Shields, Opt. Express 16, 18790 (2008). 
[74] T. Schmitt-Manderbach, H. Weier, M. Fiirst, R. Ursin, 

F. Tiefenbacher, T. Scheidl, J. Perdigues, Z. Sodnik, C. 

Kurtsiefer, J. G. Rarity, A. Zeilinger, and H. Weinfurter, 

Phys. Rev. Lett. 98, 010504 (2007). 
[75] R. J. Hughes, W. T. Buttler, P. G. Kwiat, S. K. Lamore- 

aux, G. L. Mokgarn, J. E. Nordholt, and C. G. Peterson, 

IEEE Aerospace Conference, 1803 1, 191 (2000). 
[76] J. G. Rarity, P. R. Tapster, P. M. Gorman, and P. Knight, 

New. J. Phys. 4, 82 (2002). 
[77] E. Miao, Z. Han, S. Gong, T. Zhang, D. Diao, and G. 

Guo, New J. Phys. 7, 215 (2005). 
[78] C. Bonato, A. Tomaello, V. D. Deppo, G. Naletto, and 

P. Villoresi, New J. Phys. 11, 045017 (2009). 
[79] P. Villoresi, T. Jennewein, F. Tamburini, M. Aspelmeyer, 

C. Bonato, R. Ursin, C. Pernechele, V. Luceri, G. Bianco, 

A. Zeilinger, and C. Barbieri, New J. Phys. 10, 033038 

(2008) . ^ 

[80] Perkin Elmer, http: / /optoelectronics. perkinelmer.com' 
[81] K. M. Rosfjord, J. K. W. Yang, E. A. Dauler, A. J. Ker- 

man, V. Anant, B. M. Voronov, G. N. Gol'tsman, and 
K. K. Berggren, Opt. Express 14, 527 (2006). 

[82] D. Rosenberg, S. W. Nam, P. A. Hiskett, C. G. Peterson, 
R. J. Hughes, and J. E. Nordholt, Appl. Phys. Lett. 88, 
021108 (2006). 

[83] M. Lipson, Nanotech. 15, S622 (2004). 

[84] S. L. Valley, Handbook of Geophysics and Space Environ- 
ments (McGraw-Hill Book Co., New York, 1965). 

[85] S. Passman and L. Larmore, Atmospheric Transmission 



12 



(The Rand Corporation, Santa Monica, 1956). 
[86] A. Berk, L. S. Bernstein, and D. C. Roberston, 
Tech. Rep. GL-TR-89-0122 (Air Force Geophys. Lab., 
Hanscom Air Force Base, Mass., 1989). 



[87] R. L. Fante, Proc. IEEE 63, 1669 (1975). 
[88] L. C. Andrews, R. L. Philhps, and P. T. Yu, Appl. Opt. 
34, 7742 (1995). 



